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Abstract 



An approach is developed, combining the ideas of quantum statisti- 
cal mechanics and multichannel theory of scattering, for treating statisti- 
cal systems whose constituents can possess different bound states realized 
as compact clusters. The main principles for constructing multichannel 
cluster Hamiltonians are formulated: principle of statistical correctness, 
principle of cluster coexistence, and principle of potential scaling. The im- 
portance of the principle of statistical correctness is emphasized by showing 
that when it does not hold the behaviour of thermodynamic functions be- 
comes essentially distorted. And moreover, unphysical instabilities can 
appear. The ideas are carefully illustrated by a statistical model of hot 
nuclear matter. 
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1 Introduction 



If the constituents of a statistical system strongly interact with each other, they 
can form different bound states exhibiting themselves as compact clusters. Ubiq- 
uitous chemical elements are examples of such clusters [1]. Another well known 
example is fog consisting of droplets made of water molecules. A variety of other 
examples are listed in ref. [2]. 

Despite their so widespread prevalence there is no good theory describing 
statistical clustering systems. It is possible to separate three problems in this 
description: One is the consideration of the dynamics of growing or evaporating 
clusters. This process can be portrayed by complicated partial differential equa- 
tions, similar to those that treat the dynamics of nuclei of a new phase or the 
growth and dissolution of macrodefects in nonequilibrium systems [3-5]. Another 
problem is the description of the properties of a single cluster inside equilib- 
rium matter. This can be exemplified, for instance, by the model of a cluster in 
thermostate [2]. A fullerene molecule can also be treated as a large cluster in 
thermodynamically equilibrium surrounding [6,7]. A nucleon in nuclear matter 
can be considered as a quark bag or quark cluster [8,9]. The third problem is 
to give a statistical description of a quantum system whose particles can form 
various bound states, so that in the same system different types of clusters arize. 

In this paper we address the third problem of the mentioned three, that is, the 
statistical description of a quantum system with multiple clusters. Theoretical 
treatment of such a system is the least understood and developed, as copmared to 
the first two problems, There are several objective difficulties in treating clustering 
matter. The main of them are as follows. 

Suppose that particles forming a statistical system have various bound states. 
Each bound state can be described by a Bethe-Salpeter-type equation for a z - 
particle Green function, where z is the compositeness of the bound state. This 
type of equations, even for a two-particle Green function, is very difficult to deal 
with. When there are many bound states interacting with each other as well 
as with unbound particles, then one has to deal with a large system of many 
coupled Bothe-Salpeter equations for two-, three-, and so on up to the highest 
order z -particle Green functions, plus the Dyson-type equations for single- 
particle Green functions [10]. Even to deal with a two-particle Green function is 
not as easy. And imagine that one has to operate with a ten-particle Green func- 
tion, when there occur ten-particle bound states. It is evident that in the case 
of multiple bound states the standard approach becomes practically unsolvable 
and useless. Then one simplifies the problem by inventing effective thermody- 
namic potentials. However, because of ambiguity in constructing such potentials, 
the thermodynamic behaviour of the considered system may have quite different 
features depending on the assumptions used. Morover, the fabricated thermody- 
namic potential may have no relation not only to the considered clustering matter 
but even to any statistical system. Thus, the first question that arises is: What 



3 



is a general criterion which any effective thermodynamical potential must satisfy 
for correctly representing a statistical system? 

In order to simplify the problem one often postulates that solely one type of 
bound states can exist under the given thermodynamical parameters, say temper- 
ature and density. An ensemble of given one-type bound states then is equivalent 
to a thermodynamic phase. And a transformation of one phase into another goes 
through a genuine phase transition. Evidently, such a simplification is not merely 
too rough but contradicts the general properties of a quantum system in which 
bound and unbound states can be formed simultaneously. The latter is sinony- 
mous to requiring that different types of clusters could coexist. The necessity of 
allowing for such a cluster coexistence is advocated and stressed in the present 
paper. 

As soon as we accept the possibility of cluster coexistence, another problem, 
looking unsolvable, emerges. When the particles constituting a system can form 
several types of bound-state clusters, then it is necessary to define many inter- 
action potentials between these clusters as well as between the clusters and the 
constituent particles. If the particles can make, say, ten types of clusters, then we 
need to know Cf + 10 = 55 different interaction potentials. Where the latter 
can be taken from? 

The aim of this paper is twofold: To formulate the general principles of a 
correct statistical description of clustering matter and to illustrate these principles 
by examples. 

2 Principle of Statistical Correctness 

Wishing to analyse thermodynamic properties of a statistical system, we need to 
define a thermodynamic potential. Writing the latter, one usually invokes some 
simplifications leading to an effective thermodynamic potential Q((p) containing 
a set tp = {tpi} of auxiliary functions depending on space and/or thermodynamic 
variables. What necessary to keep in mind, first of all, is that not each effective 
potential can have sense, however reasonable it may look. A thermodynamic 
potential, to be accepted as such, must satisfy the properties formulated below. 

Property 1. Statistical Representability. 

An effective thermodynamic potential fl(tp) represents an equilibrium sta- 
tistical system if and only if it has the Gibbs form 

f%) = il[H{<p)] (1) 

defined as 

n[H] = -TlnTr e~^ H , (2) 

where f3 = T~ x , and the dependence on auxiliary functions comes only through 
an effective Hamiltonian H(ip) . Such a thermodynamic potential is called 
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statistically representable. 

Thus, if one invents an effective thermodynamic potential, even pronounc- 
ing seemingly plausible words, this does not necessary mean that the invented 
potential describes some statistical system. If this potential is not statistically 
representable it corresponds to no equilibrium statistical system. For example, a 
thermodynamic potential in the excluded-volume approximation is not statisti- 
cally representable. 

Property 2. Thermodynamic Equivalence. 

A statistical system defined by a Hamiltonian H is thermodynamically 
equivalent to a system modeled by an effective Hamiltonian H((p) if and only 
if their thermodynamic potentials Q and £l((p) are statistically representable, 
i.e. 

n = q[h], n((p) = n[H((p)], (3) 

and are equal to each other, 

n[H] = n[H(<p)]. (4) 

The corresponding Hamiltonians are called thermodynamically equivalent. 

For the case of infinite matter, condition (4) can be weakened by requiring 
the validity of the asymptotic, in the thermodynamic limit, equality 

Jim 1 (Q[H\ - Q[H(<p)]) = 0. 
Property 3. Statistical Correctness. 

The necessary condition for an effective Hamiltonian H(tp) to be thermody- 
namically equivalent to a Hamiltonian H , not depending on auxiliary functions 
ip , is the variational equality 

= 0, (5) 

where the variation over ip implies the set of variations with respect to each tpi , 
and the average of an operator A is 

~ Triexp{-/?i%)} 
[ ' ~ Trexp{-(3H(<p)} ' 1 ' 

Equality (5) is an evident consequence of the property of thermodynamic equiv- 
alence. We shall call (5) the condition of statistical correctness. 
Basing on these properties, we formulate the 

Principle of Statistical Correctness: 
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An effective thermodynamical potential is statistically correct if it is statisti- 
cally representable with an effective Hamiltonian satisfying the condition of sta- 
tistical correctness. 

Suppose that the thermodynamic potential Q = fl(tp) is statistically repre- 
sentable with the Hamiltonian 

H(<p)=E (7) 

i 

in which E is the energy operator, /ij is a chemical potential, and Ni is 
a number-of-particle operator. If this thermodynamic potential is statistically 
correct, then condition (5) guarantees the validity of the thermodynamic relations 

P ~~dV ~ ~V' 
Op \ 



dp 

S= df 




J2ViPi), (8) 



for the pressure p , energy density e , entropy density s , and particle density 
Pi ; the volume of the system being V . 

For the effective thermodynamic potential not satisfying condition (5) the 
thermodynamic relations (8) would be broken, that is, the values of the quantities 
on the left-hand side of (8) calculated as the corresponding derivatives or as 
statistical averages would be different. Such an inconsistency in (8) would signify 
that one should not trust to predictions derived from an effective thermodynamic 
potential which is statistically incorrect. 

The importance of sustaining the self-consistency in the thermodynamic rela- 
tions (8) for studying the thermodynamics of effective models was emphasized by 
Zimanyi et al. [11]. Requiring the validity of (8) for an effective thermodynamic 
potential containing unknown correcting functions yields a complicated system of 
nonlinear differential equations in partial derivatives with respect to the variables 
T, V , and /ij . This system of equations is to be complimented by boundary 
conditions. Such a system has no unique solution. To extract somehow the lat- 
ter, one needs several additional heuristic assumptions and fitting parameters. 
Contrary to this, the principle of statistical correctness formulated above is, as 
we shall show in the next section, much simpler to deal with and gives a unique 
solution. 

Also, the principle of statistical correctness is more general leading to (8) 
but not conversely. Moreover, if one finds correcting functions directly from the 
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first-order relations (8), this does not guarantee the validity of the second order 
relations for the specific heat 

Cv = % = £ (<^ 2 > - <^> 2 ) 0) 
and the isothermic compressibility 

^-H^r -<*>')■ « io) 

At the same time, the condition of statistical correctness (5) insures that (8), (9), (10) 
and other analogous relations of arbitrary order hold true. 



3 Principle of Cluster Coexistence 

If the particles of a quantum system form different bound states, then the lat- 
ter can influence thermodynamic properties. To take this influence into account, 
one has to allow for the existence of these bound states, estimating their rela- 
tive contribution to the properties of the system. This is equivalent to saying 
that it is necessary to allow for the coexistence of different clusters, calculating 
their statistical weights in the analyzed properties. The necessity of taking into 
account both unbound and various bound states seems so natural from the gen- 
eral quantum-mechanical point of view that it could look excessive repeating it, 
if it would not be so common meeting in literature statistical models in which 
different cluster states are prohibited to coexist, being treated as pertaining to 
different thermodymanic phases. Because of this frequent confusion of quantum 
states with statistical states we feel it is worth clarifying this question once more. 
A quantum bound state corresponds to a cluster but not to a thermodynamic 
phase. Although it may happen that under some conditions, say at low tempera- 
ture, one type of clusters prevails while at high temperature another type becomes 
dominant, this does not mean that there is no admixture of other types of clus- 
ters among those of a predominant type. Predominance of one type of clusters 
is not the same as the complete prohibition to exist for other types. What types 
of clusters and in which proportions can coexist in a particular thermodynamic 
phase is, figuratively speaking, to be decided by the system itself, which implies 
that each system tends to a state of maximal thermodynamic stability, and the 
concentrations of clusters are to be defined by stability conditions. The latter can 
be essentially spoiled by the prohibition for clusters to coexist and can lead to 
incorrect thermodynamic behaviour, for instance, to the appearance of spurious 
phase transitions or to the change of order of genuine phase transitions. 

The most general and, probably, apparent way of understanding the structure 
of a Hamiltonian corresponding to a system of coexisting clusters is from the 
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point of view of the multichannel theory of scattering [12] considering each type 
of bound states, i.e., of clusters, as a reaction channel of interacting particles. 

The methods of the abstract multichannel theory of scattering may be applied 
to physical systems of different nature, in which the constituent particles can form 
various bound states. Below we briefly delineate the main general ideas that could 
be used for any kind of such systems. 

Let a system be defined by a Hamiltonian H which is a selfadjoint operator 
acting in the Hilbert space H , called the space of quantum states. Assume that 
the constituent particles of the system can form bound states. Thus, in a quark- 
gluon system various hadron states can be formed. Enumerate all possible types 
of bound states by the index % , where % — 1 stands for unbound particles. Each 
type of bound states is individualized by a set of corresponding characteristics, 
such as the compositeness number Zi showing the number of bound particles, 
effective mass , and a set of quantum numbers like spin, isospin, colour, 
baryon number, and so on. All quantum states associated with the same type 
of bound states, indexed by i , compose the subspace Hi C H . In other 
words, Tii = P{H is a projection of H . Generally, different quantum states can 
be made orthogonal to each other. This means that the subspaces Hi can be 
considered as mutually orthogonal, that is, for the reducing projections p one 
has PiPj = 5ij . Such pairwise orthogonal projections p are called the channel 
projections, and the related subspaces Hi are termed the channels. A quantum 
state pertaining to the channel Hi is an % -channel state. 

The set of channels Hi is complete if all possible % -channel states span the 
whole space of quantum states H . Then the latter is written as the direct sum 
H = ®{Hi . This is equivalent to the resolution of unity 1 = Si P • 

Let the time evolution of quantum states in a channel Hi be defined by a 
selfadjoint operator Hi . This implies that the channel Hi is invariant under 
the action of Hi . Because of the pairwise orthogonality of the channels, the 
operators Hi are pairwise commuting, [Hi, Hj] = . Each of such operators 
Hi is called the channel Hamiltonian. The channel system is a system {Hi} 
of the channel Hamiltonians, each of which acts in the corresponding channel 
Hi . The sum J2i Hi + const 1 of the channel Hamiltonians, attributed to the 
channel system, may be named the multichannel Hamiltonian. The latter, by 
construction, acts in the space of quantum states @{Hi . 

In this way, the description of an ensemble of particles forming bound states 
can be done by constructing the corresponding channel system. To this end, one 
has to classify different types of states and to define the related channel Hamil- 
tonians Hi . The sum J2i Hi + const 1 gives the Hamiltonian of the channel 
system. It is worth noting that the multichannel Hamiltonian does not necessarily 
coincide with the exact Hamiltonian H , but rather gives a physically reasonable 
first approximation. In particular calculations, one can, if necessary, resort to per- 
turbation theory starting from the multichannel Hamiltonian. Though in many 
applications already the channel approximation yields quite accurate results. 
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Summarizing what is said above we formulate the 
Principle of Cluster Coexistence: 

If the particles of an equilibrium statistical system can form clusters of differ- 
ent types, then all such cluster types can coexist with the probabilities defined by 
the condition of thermodynamic stability. 

Let us concretize how one can define such channel probabilities. The multi- 
channel Hamiltonian, in general, reads 

H = Y,H l + CV ) (11) 

i 

where Hi is an % -channel Hamiltonian and CV , a nonoperator term. 
Take the channel hamiltonians, Hi , in the mean-field approximation 

H l = Y J <k)a\(k)^(k) (12) 

k 

with an effective spectrum 

Ui(k) = Ki(k) + Ui- fii, 

in which k is the absolute value of the momentum k ; Ki(k) is the kinetic- 
energy term; Ui , a mean field; and \ii , the chemical potential of the % -type 
cluster. The field operator a, (A;) is a column in the space of internal degrees 
of freedom, such as spin, flavor, and so on. For the mean-field Hamiltonian (12) 
the momentum distribution 

(a\(k)ai(k)) = Oni(k) 

is easily calculated; Q being a degeneracy factor, and 

m(k) = [exp{Pui(k) =F I}]" 1 

is the Bose- or Fermi function depending on the upper or lower sign, respectively. 
The average density of the % -type clusters is 



rii{k)d k ■ (13) 



Note that the isotropicity of the system is assumed here, because of which rii(k) 
depends on k = | k \ ■ Consequently, in (13) we could write 



rii(k)d k— 47r J ni{k)k 2 dk; 
but for the sake of brevity we prefer the former notation. 
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The density of % -clusters (13) factored by the compositeness number zi 
gives the density of particles ZiPi that are bound in the clusters of the % -type. 
The total average density of particles is 

P = H z iPi- ( 14 ) 

i 

The statistical weight of each channel is characterized by the channel probability 

id = Zi — , (15) 
P 

of the i -type clusters. By definition, eq. (15) enjoys the conditions 

0<Wi<l, £>i = l. (16) 

i 

The compositeness number of unbound particles, Zi , is assumed to equal one. 

The chemical potentials of i -clusters, entering into (15), can be expressed 
through the given thermodynamic variables with the use of the equilibrium con- 
ditions considered. Thus, if the average density of particles (14) is given, then 
the i -cluster chemical potentials pi are related to that of unbound particles, 
p , by the equilibrium condition 

— = p (p = const). (17) 

In the case when the average baryon density 

n B = J2BiPi (18) 

i 

is conserved, where Bi is the baryon number of a cluster in an i -channel, then 
the equilibrium condition reads 

^ = p B {n B = const), (19) 

where p B is the baryon potential of matter. 

The cluster probabilities in (15), together with the kind of equilibrium consid- 
ered, are defined as functions of thermodynamic variables. Finally, it is necessary 
to check that the considered equilibrium is stable, so that the stability conditions 

C v > 0, k t > (20) 

for the specific heat (9) and isothermic compressibility (10) are fulfilled. 

In this way, all cluster probabilities in (15) can be calculated self-consistently, 
thus showing the proportions in which various types of clusters are intermixed in 
matter. 
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4 Principle of Potential Scaling 



In order to complete the definition of a system containing different types of clus- 
ters, we need to specify the interaction potentials $ij(r) of these clusters. If 
there is a number of such cluster types, then we have to define numerous poten- 
tials &ij{r) . For example, as is mentioned in Introduction, in the case of 10 
cluster types we need to have 55 such potentials. And for 100 cluster types 
we would need already 5050 interaction potentials. Where could we take them 
from? 

The idea of how the interaction potentials could be connected with each other 
comes from the form of the equilibrium conditions (17) and (19). By analogy to 
these conditions we may conjecture that the interaction potentials are scaled, by 
means of the corresponding compositeness numbers, to a universal function as 
follows: 

= $(r). (21) 

ZiZj 

Such a relation, as is clear, can be sensible only if the interaction potentials are 
of similar nature. To be more precise, we shall say that the interaction potentials 
are in the same universality class if and only is they can be scaled to one universal 
function, as in (21). As a counterexample we may adduce a pair of potentials 
one of which decreases, as r — > oo , and another increases. Certainly, these 
potentials cannot be in the same universality class. Only the clusters of similar 
nature possess the interaction potentials that can pertain to one such a class. 
Combining what is said above, we come to the 

Principle of Potential Scaling: 

The clusters of similar nature interact with each other through potentials from 
the universality class, so that these potentials are scaled by the corresponding 
compositeness numbers to a universal function. 

The scaling relation (21), defining a universality class, can be derived from 
the following reasoning. Consider the channel reaction 

m + n + j^i+j (22) 

with the corresponding equation for the compositeness numbers 

z m -\- z n -\- Zj — Zi -\- Zj. (23) 

The reaction (22), with eq. (23), signifies that the m -cluster and n -cluster fuse 
together in the presence of a j -cluster. Assume the additivity of interactions 
with respect to the j -channel, that is, 

* mj (r) + * nj (r) = * iS (r). (24) 
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Then, from (23) and (24) it is straightforward to derive (21). Therefore, the 
scaling relation (21) can be interpreted as the result of the additivity of interac- 
tions (24) under condition (23). Making this statement more general, but less 
accurate, we may say that the scaling relation (21) is a manifestation of two 
conservation laws: conservation of energy and conservation of particle number 
during the reactions of fusion and decay. 

Thus, if at least one of the interaction potentials is known, say $ mn (r) , then 
by means of the scaling relation (21), all others can be expressed through this 
reference potentials as 



This solves the problem of multiple interaction potentials making the multichan- 
nel approach to clustering matter completely defined. 



To illustrate the approach, we need to apply it to a particular clustering substance, 
it is possible to distinguish two quite different cases: One, when the consideration 
can be limited by several cluster types, say, by about ten of them; and another, 
when there exists a multitude of various cluster types. The latter case has to do 
with substances like fog which consists of droplets whose compositeness numbers 
range from one to many billions. A good example of the former case is nuclear 
matter which can be characterized by a limited number of clusters. In this way, fog 
is a more complicated system, which we shall consider in a separate publication. 
And here we analyse the more simple case of hot and dense nuclear matter, 
paying the main attention to the process of evaporation and condensation of 
hadron clusters, that is, to the deconfinement-confinement transition. 

The necessity of resorting to statistical models for nuclear quark-hadron mat- 
ter is due to the fact that perturbative quantum chromodynamics does not pro- 
vide information on the equation of state in the whole region of thermodynamic 
parameters. Especially little can be said about the most interesting region of 
deconfinement-confinement transition. There exist lattice calculations which, 
however, are reliable only for the case of zero baryon density, but not for its finite 
values [13]. 

The overwhelming majiority of statistical models for deconfinement have been 
based on the assumption that the latter is a phase transition between two pure 
phases, the clustered hadron phase and the phase of unbound quarks and glu- 
ons called the quark-gluon plasma. But this assumption contradicts to many 
computer simulations and also to some analytical estimates. For example, in- 
tensive numerical simulations for lattice quark-gluon plasma [14] and pure gluon 
plasma [15] revealed nontrivial effects due to strong particle correlations at tem- 
peratures above the deconfinement temperature T d , which has been interpreted 
as the existence of hadronic modes, even at T > T d . The deconfinement was 




(25) 
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found to be a rapid crossover, but not a genuine phase transition [14]. The high- 
temperature phase contains fluctuations being colour singlet modes, hadronic in 
character. This is because the poles and cuts in the linear responce functions of 
the hadronic phase go over smoothly into those of the high-temperature one, so 
that all low-temperature hadrons have analogs in the high-temperature phase. 
The low-and high-temperature phases both have confining characteristics, but 
the effect of confinement upon thermodynamic properties becomes less and less 
significant as temperature increases. 

The confinement-deconfinement transition can be compared with the insulator- 
metal transition or with the ionization, since hadrons are nothing but the bound 
states of quarks and gluons, or small droplets of quark-gluon plasma [16-18]. In 
insulating solids, below the transition point to a metal, the conductivity is not 
strictly zero, since the ionization energy is finite and ionization can locally provide 
some free electrons. Similarly, in the hadron matter below the deconfinement tem- 
perature quarks can be separated. Both the conventional Mott transition in solids 
and the deconfinement transition in hadronic matter thus lead from a regime, in 
which the binding can locally be broken by ionization, to one where it is globally 
removed by a collective screening of the binding force. 

The spatial structure of correlation functions, obtained in lattice numerical 
simulations for hot quark-gluon plasma, is very similar to the structure of the cor- 
responding zero-temperature functions [14]. The thermodynamic characteristics, 
such as pressure and energy density, aslo display strong nonperturbative effects 
persisting till about 2Tj , as is found in the lattice simulations and discussed 
in recent surveys [19,20]. The fact that some lattice simulations with dynamical 
quarks displayed jumps of thermodynamic characteristics at Td can be explained 
[21] as merely due to the finite size of the lattice, since the transition becomes 
less and less abrupt as the lattice size inreases. Lengthy runs show no evidence 
for metastability thus suggesting that there is no sharp transition, but only a 
crossover phenomenon. 

The coexistence scenario, based on numerical simulations, is supported also 
by some analytical calculations. For instance, the density-density correlation 
functions in the Nambu-Jona-Lasinio model contain poles and cuts that are the 
same at all temperatures [22]. Therefore meson modes do exist above as well as 
below transition temperature. Analogously, quarks and gluons should also exist 
in the low-temperature as well as in high-temperature phase. The graduate 
change in the excitation spectrum from hadronic states to quarks and gluons 
and the survival of hadronic modes in the high-temperature phase appear also 
in the magnetic-current approximation [23] and in the instanton-liquid approach 
[24]. This suggests the following picture of the quark-hadron coexistence. The 
correlations between quarks persist above the deconfinement temperature forcing 
some of them to correlate into colour singlets. As the quarks are moving in 
the heat bath, the strings connecting them for colour neutrality are constantly 
breaking and reforming, which can be interpreted as hadrons going in and out 
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of the heat bath. This picture has a formal resemblance to the string-flip model 
[25] although with time-like strings. 

Thus, the concept of cluster coexistence is absolutely natural from the point of 
veiw of the multichannel scattering theory and is supported by lattice numerical 
calculations. The cluster coexistence is somewhat similar, although not identical 
to heterophase coexistence common for many statistical systems, both these types 
of coexistence leading to precursor, or pretransitional, phenomena [26]. 

There is a variety of different statistical models trying to describe deconfine- 
ment in nuclear matter (see reviews [27,28]). The majority of these models, with 
a few exceptions (e.g. [29-31]), do not take into account the possible coexistence 
of clusters. This is why the predictions of such models are in disagreement with 
lattice numerical results. 

Consider the clustering nuclear matter whose elementary constituents are 
quarks, antiquarks and gluons. These can form bound states corresponding to 
hadron clusters. The total set {i} of the indices enumerating the channels can 
be separated into two groups, and {i} z . The first group is related 

to unbound constituents, quarks, antiquarks, and gluons whose compositeness 
number Zi — 1 . The second group corresponds to hadron channels representing 
bound states with compositeness number Zi > 2 . Respectively, the density (14) 
may be written as the sum 

p = pi+p z , p! = Y,Pi, Pz = ^ZiP h 

in which p\ is the density of unbound particles and p z is the density of particles 
in bound states. 

Accept the effective Hamiltonian (11) with the channel Hamiltonians given 
by (12). To define the effective spectrum LUi(k) in (12), we have to concretize 
the mean field Ui . The quark-gluon plasma mean field U\ can be defined as 

Z7i = U(p) = pj V(r)s(r)d r, (26) 

where V(r) is a confining potential and s(r) , a screening function. The ne- 
cessity of including a screening correlation function into the mean-field approx- 
imation is the general requirement for any statistical system with nonintegrable 
interaction potentials [32]. The approximation (26) is the correlated Hartree ap- 
proximation. The confining potential is not integrable since its mostly often used 
representation has the power-law behaviour 

V(r) = Ar u (0<v< 2). (27) 

The screening function 
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can be scaled by the mean interparticle distance a . Then, with the notation 

roc 

J 1+u = An A / s(x)x 2+u dx (28) 
Jo 

for the effective intensity of interactions J measured in energy units, the mean 
field (26) becomes 

U(p) = J l+U p-^ 3 . (29) 
The plasma mean field (29) has the following asymptotic properties: 

oo, p — > 0, 

(30) 

0, p — > oo. 

The upper line here tells that quarks and gluons cannot exist as free particles 
outside dense nuclear matter - this is what is called the colour confinement. 
The lower line in (30) shows that the interparticle interaction decreases with 
the decrease of the average distance between particles - this is the so called 
phenomenon of asymptotic freedom. 

The mean field for an i -channel corresponding to a bound hadron state can 
be written as 

U i ='£* ij p j + z i [U(p)-U(p x )], (31) 
{jh 

where the summation is over the bound states, and p z is the density of particles 
in bound states. The first term in (31) desribes the interaction of a cluster in an 
% -channel with all other hadrons; the mean interaction potential being 

= J ^( r ) d r, %(r) = Vij{r)8ij{r), 

where Vij(r) is a bare interaction potentials between the clusters of the i - 
and j -types, and s^(r) is a screening correlation function, so that $ij(r) is 
called the screened, or effective, potential. The second term in (31), in the square 
brackets, models the interaction of a cluster of the % -channel with unbound 
plasma states. 

Invoking the principle of potential scaling (21) for the effective potentials 
&ij(r) , we may write 

= Zizfi, $ = J $(r)d r, (32) 
where $(r) is some reference function. Then (31) reduces to 

Ui = Zi$p, + z,J 1+v [p-^ - p^ 3 ) . (33) 

At this point it is necessary to stress that the interaction between unbound 
particles and these between hadron clusters are of different nature, the former 



U(p) -> 



15 



growing while the latter diminishing with increasing interparticle distance. This is 
equivalent to saying that they are from different universality classes, thus, cannot 
be scaled to the same reference function. Therefore, the reference function $(r) 
in (32) cannot be interpreted as an interaction potential for unbound plasma 
constituents. Consequently, with these two universality classes, we have to keep 
two parameters, $ and J , which are to be chosen independently. 

The Hamiltonian (11) must satisfy the principle of statistical correctness (see 
Sec. 2). The role of auxiliary functions in (11) is played by the densities p 
and p z , which are functions of temperature and baryon density. Therefore, the 
condition of statistical correctness (5) reads 



0, — = 0. (34) 





Note that instead of p and p z , as auxiliary functions, we could take p and 
pi or pi and p z . The condition (34) does not depend on this choice, since 
p = pi + p z . Substituting (11) and (12) into (34) gives 

i Pl Sp bp ' Y %s Pz s Pz 
Taking account of (31) and (33) yields 

The solution of these equations is straightforward and, up to a constant, it is 

C = ^-J^» (pWs _ p W3) _ U pl (35) 

In this way, the multichannel Hamiltonian (11) is completely defined. Em- 
phasize that the term CV cannot be omitted since its presence provides the 
validity of the principal of statistical correctness. Only retaining this term makes 
it possible to find the correct thermodynamic behaviour of the clustering matter. 

With the Hamiltonian (11), the pressure is 

P=Eft» p i = ±Tj^J\n[l±n i {k)]dk, (36) 
and for the energy density one has 

(if ~* 

e = y^e i) £i = —1— J [^k) + p} ni {k)d k ■ (37) 

Taking the kinetic-energy term in the relativistic form Ki(k) = \Jk 2 + m, 
where m 8 is the corresponding mass, we get the spectrum 



2 



LUi(k) = ^Jk 2 + ml + Ui- in. (38) 
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Before analysing in detail the thermodynamic behaviour of the system, let 
us notice that it has the following properties. In the case when T — > and 
tib — > , all particles are condensed into hadron clusters. And when T — > oo , 
at any fixed n# , only unbound states survive. Consider the high-temperature 
case more accurately, since in this limit we may compare the results with the 
available perturbative calculations in quantum chromodynamics. 

When T — > oo , then only quarks, antiquarks, and gluons remain, so that the 
set of unbound states {i}i is {q, q, g} . The pressure (36) and energy density 
(37) become 

p~Y,Ph e~Y, £ i- ( 39 ) 

At temperatures much higher than the quark masses the latter can be neglected, 
because of which in (36) and (37) we may substitute 

rii(k) — > [cxp{(3(k - /ij)}=F] _1 , (40) 

the chemical potentials being 

Hq = -Hq =H, Vg = 0. (41) 

Thus we come to 

d f°° k 3 dk _ 

Pt ~6^Jo exp^A;-/^)}^' El ~ Pl ' 1 } 

Take into account that the degeneracy factors of quarks and antiquarks are the 
same, ( g = (g , and that these particles are fermions while gluons are bosons. 
Then an exact integration yields 

From here, for the pressure in (39) we have 



^s(i^) r4+ !" 2r2 ( i+ 2^i- < 43 » 



For the specific heat we find 



C v - 



15 W q + ^ 9 J + 2(3/i 2 + 7r 2 T 2 ) C ' 9 ' 1 } 



The pressure (43) is to be compared with that obtained by perturbation theory 
in quantum chromodynamics (see Appendix). Due to the relations 

( q = 2xN f xN c , ( g = 2x (N 2 C - 1) (45) 
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for the degeneracy factors of quarks and gluons, where Nf and N c are the 
numbers of flavours and colours, respectively, we make it sure that the high- 
temperature pressure (43) asymptotically coincides with the perturbative pressure 
in QCD , as t — > oo . 

Expression (43) can be simplfied further if we substitute there the chemical 
potential as a function of T and n B , which can be found from the formula for 
the baryon density 

n B ^ \{p q - p 9 ) = / [n q (k) - n q {k)]d k . (46) 

With eq.(40), this gives 

n B ~^(v? + «>T°). (47) 
From (47), we get the chemical potential 



(T-.oo). (48) 



For the density of quarks and gluons we find 

Pq c |%3)T 3 , Pg ~ %(3)T 3 , 

4:71 7T Z 

where ((3) = 1.20206 . Using (48) and introducing the notation 

C = \c q + Cg, (49) 

we come to the conclusion that the high-temperature behaviour of the system 
asymptotically reduces to that of the Stephan-Boltzmann plasma, 

P - Psb, e ~ e SB , C v ~ C SB , 

for which 

7T 2 7T 2 27T 2 

Psb = -(T\ £ sb = -CT\ C SB = —CT\ (50) 

with the factor (49). In what follows it will be convenient to present the results 
in a dimensionless form scaling them by means of the Stephan-Boltzmann ex- 
pressions in (50). For the number of colours N c = 3 , we have ( q = 6Nf and 
( g = 16 . Therefore, the factor (49) is 

21 



C = 16(l + S AT / ), 



which is to be substituted into (50). Then, e.g., the pressure becomes 

7T 2 / 21 



^ = Io-v 1 + 32^) T 
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6 Analysis of Thermodynamic Characteristics 



To accomplish explicit calculations of thermodynamic characteristics, we need to 
specify three parameters, J, $ , and v . The first of them, J , describes 
the intensity of nonperturbative interactions in the quark-gluon plasma. We can 
evaluate J in several ways. The simplest way is to remember that nonper- 
turbative effects in the quark-gluon plasma are commonly associated with the 
so-called bag constant, B , for which one uses different values so that B 1 ^ lies 
in the interval between about 150 to 300 MeV . We can accept for J the 
value from the middle of this interval, that is J = 225 MeV . 

The constant <£> , according to (32), can be chosen as $ = $33/9 with 
$33 given by $33 = / Vs3(r)s3s(r)dr , where V^r) is the interaction potential 
between the three-quark bound states, that is between nucleons. The nucleon 
interaction potentials are well known from scattering experiments. There are 
several representations for these potentials. For nonintegrable hard-core poten- 
tials one has to take into account the correlation function S33(r) , while for 
integrable soft potentials 533(7") m 1 . Among many known nucleon-nucleon 
potentials, we prefer the Bonn potential [33], which provides an accurate descrip- 
tion of nucleon scattering and has sufficiently simple analytic form. Following 
the common consensus that thermodynamics of nuclear matter does not depend 
on the mutual orientation of spins of scattering nucleons, we average over spin 
directions, and assume also that the interaction between any pair of nucleons, 
whether these are protons or neutrons, is the same. The so-called cut-off terms 
of the Bonn potential can be neglected since they start playing an essential role 
only for very short distances < 0.1 fm , which would correspond to the baryon 
density n B > 10 3 n _B at which it would be hard to expect that any nucleons 
could survive. Here and in what follows n s = 0.167 fm~ 3 is the normal baryon 
density. The parameter $ calculated in this way is $ = 35 MeV fm 3 . 

The constant v characterizes the power-law behaviour of the confining po- 
tential (27). The most frequently accepted cases are those of linear and harmonic 
confinement. Actually, our results do not change qualitatively for v between 
these two possibilities. Let us consider, for example, the harmonic confinement. 
As we have checked, our results practically do not change quantitatively in the 
region 1.5 < v < 2 . Thus, in what follows j/rj2. 

Finally, we have to concretize the reaction channels that will be included 
into consideration. In principle, the developed approach permits us to include 
any number of particles. We have analysed many variants which, as we found, 
demonstrate similar behaviour. Not to overload this paper, here we limit ourselves 
by the particles listed in the Table. We included into consideration multiquark 
states, although their status is not yet absolutely clear, because there has been 
an intensive discussion about their possible presence in nuclei (see reviews [34- 
36]) motivated by the Baldin interpretation of the commulative effect [37]. The 
chosen mass of the 6 -quark state corresponds to an average over the masses of 
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several light dibaryons that are claimed to be observed in experiments [38]. The 
9 - and 12 -quark parameters are elicited from the bag model calculations [39]. 

We have analysed in detail the behaviour of the main thermodynamic charac- 
teristics of the model as functions of temperature 6 = k B T and relative baryon 
density n B /n B ■ Fig.l shows the probability of the plasma channel 

Wpl = ^{pg+ Pu+ Pu+ Pd+ Pd) ■ 

The probabilities of hadron channels are displayed in the following figures: in 
Fig. 2, the pion-channel probability 

2 / 

W n = - {p n + + p w - + p n o) , 
P 

in Fig. 3, the summarized, excluding pions, probability of other meson channels 

2 / 

W V pw = -{P V + Pp+ + Pp- + Pp" + Plu) , 

in Fig. 4, the nucleon-channel probability 

W 3 = - {Pn + Pn + Pp + Pp) , 

in Fig. 5, the probability of the six-quark channel 

6 < 

W 6 = - (pQq + p 6 g) , 

and in Fig. 6, the probability of the six-quarks in the Bose-condensed state. The 
probabilities of heavier multiquark channels are not depicted since they are neg- 
ligibly small, 

9 12 

Wg = - (p 9 q + P9q) < 10~ 3 , Wl2 = — (Pl2q + P\2q) < 10~ 5 . 

P P 

Figs. 7-9 present some channel probabilities at zero temperature for a larger region 
of the relative baryon density. As is seen, the channel probabilities are continuous 
functions of their variables. The location of the deconfinement transition can be 
conventionally associated with the maxima of some of the channel probabilities. 
Fig. 10 demonstrates the pressure 

P = ^Pi-C ) Pi = ±T-^^ J In [l±7ii( k)] d k ; 

and Fig. 11, the energy density 

e = Y^Ei + C, Ei = J u)i{k)rii{k)d k +Bip B Pi- 
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These functions are quite monotonous but their ratio 



V --c 2 

having the meaning of the effective sound velocity squared, exhibits in Fig. 12 a 
maximum around = 160 MeV . The temperature dependence of c 2 e ^ at 
jig = agrees with that reconstructed from lattice data. 

The specific heat in Fig. 13 is also monotonous, but the reduced specific heat 

T de 

av = 7df 

in Fig. 14 again displays the maximum that can be associated with the location of 
the deconfinement line. The similar maximum exists for the compression modulus 

-i dp 
k t = n B - — 
dn B 

in Fig. 15. 

These thermodynamic characteristics demonstrate that the transition from 
bound hadron states to unbound quark-gluon states is a gradual crossover. All 
these states are different channels of the same quantum system. Each channel is a 
subspace of the total Hilbert space of quantum states. Therefore, all channels, as 
possible states, always exist. But their statistical weights, defined by the channel 
probabilities, are different and change with varying thermodynamic parameters 
such as temperature and baryon density. The channel probabilities depend also 
on the Hamiltonian parameters, for instance, masses. Thus, if we take the mass 
of a six-quark state m 6 > 2000 MeV , then the probability of this channel 
drastically drops down. Generally, among the members of the same family of 
states only those with lower masses are the most important. The states with 
heavier masses have, as a rule, much smaller statistical weights. 

Note also that the temperature behaviour of thermodynamic characteristics 
at nonzero baryon density is modified as compared to that at zero baryon density: 
Increasing n B smoothes the transition. Therefore one has to be quite cautious 
about those predictions for the processes at nonzero baryon density, which have 
been made basing on the processes at n B = . 

In conclusion, we would like to emphasize that the main in this paper are not 
particular numerical predictions, though, we think, they are quite reasonable, 
but rather the whole picture and the ideology which the latter is based on. The 
basis of the model suggested rests upon the multichannel interpretation of differ- 
ent bound states and on the principle of statistical correctness permitting us to 
construct a correct multichannel Hamiltonian. 

Within the framework of this approach many details can be changed and some 
improvements may be made. For example, instead of the Bonn potential we could 
opt for some other effective nucleon-nucleon interactions, we could include more 
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kinds of bound states, and so on. Also, we have not discussed here the role of 
strangeness whose generation can be important as a diagnostic tool for the study 
of the quark-gluon plasma signals in realistic nuclear collisions [40] . The work on 
these problems is under progress. Changing particular details of the model might 
slightly alter some quantitative values, but the qualitative picture should remain 
the same. According to this picture, the deconfmement is a gradual crossover from 
hadron states to the quark-gluon plasma in agreement with the lattice numerical 
calculations [14]. 

7 Discussion 

In this paper we suggested a general approach for treating complex systems whose 
constituent particles can form different bound states realized as compact clusters. 
The approach is based on three main notions: principle of statistical correctness, 
principle of cluster coexistence, and principle of potential scaling. The last two 
principles are compulsory for a self-consistent treatment of clustering matter. 
The disregard of the first principle, that of statistical correctness, can sometimes 
be admissible when one is looking for a very rough approximate picture of clus- 
tering. However, if one wishes to get a really accurate description, the principle 
of statistical correctness must be satisfied. Moreover, neglecting this principle 
leads in some cases to the appearance of unphysical instabilities around phase 
transition points [28]. 

To illustrate the importance of using the princliple of statistical correctness, 
let us consider a simple clustering system consisting of two species, of unbound 
gluons and of glueballs that are bound gluon clusters. Include into consideration 
five lightest glueballs with the masses found in the bag-model calculations [41]. 
Take into account that the gluon as well as glueball chemical potentials are zero, 
/li — . Consider, for concreteness, the system with the SU(2) symmetry 
when N c = 2 and the gluon degeneracy factor is ( g = 6 . And let us calculate 
thermodynamic characteristics of such a system for two cases: when the principle 
of statistical correctness is satisfied and when it does not. In the second case we 
simply put zero the correcting term CV in the Hamiltonian (11). The results 
of calculations are presented in figs. 16-21, where the specific heat, pressure, and 
energy density are given in relative units, with respect to the corresponding quan- 
tities of the quarkless Stefan-Boltzmann plasma. The relative energy density and 
pressure are compared with the lattice numerical data, fitting the parameters of 
the model so that to have deconfmement at T d = 220 MeV . The qualitative be- 
haviour of thermodynamic characteristics is similar in both cases, in the presence 
of the correcting term CV and without it. At low temperatures, practically all 
gluons are clustered into glueballs. The glueball channel probability wg drasti- 
cally drops down at the deconfmement temperature (fig. 16). The gluon channel 
probability w g behaves oppositely to wq , as is shown in fig. 17. The deconfine- 
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merit is a second order transition with a sharp peak of specific heat (fig. 18). The 
comparison of the relative energy density with lattice numerical data for the case 
of the model without the correcting term is given in fig. 19, as compared with the 
data of ref. [42], and in fig.20, as compared with those from ref. [43]. The relative 
energy density and pressure for the model satisfying the principle of statistical 
correctness are presented in fig. 21, compared with the lattice numerical data [43]. 
As is clearly seen from these comparisons, the model without the correcting term 
does not fit well the numerical data, while the statistically correct model is in 
perfect agreement with the latter. 



Appendix. High— Temperature Limit 

In section 5 we considered the high-temperature behaviour of pressure stating 
that it coincides with the high-temperature asymptotic expression for pressure 
in quantum chromodymanics. Here we show this explicitly. 

The effective coupling parameter g = g(T) in quantum chromodynamics is 
given by the equation 

AT) 



{llN c -2N f )\n{T/A) , 

where Nf is the number of flavours and A is a scaling parameter, A « 
200 MeV . The coupling g(T) — > as T -> oo . Therefore, perturbation 
theory in powers of g becomes admissible. Perturbative expansions for the free 
energy have been obtained [44,45] to order g 5 . The fourth and fifth orders of 
these expansions contain an arbitrary renormalization scale. The QCD pressure 
can be written [17] as 

V = AT\ 
where the factor A , as T — > oo , is 

A ~ A + A 2 g 2 + A 3 g 3 , 

with the coefficients 

A ^ 

An = — 

45 
Ao = - 



iVr — 1 + -N c N f + 15N c N t 

4 * } 2vr 2 T 2 V 2vr 2 T 2 J 
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N 2 - 1 / 1 1 1 u 2 \ 3/2 

3 127T \3 6 3 2tt 2 T 2 ) ' 

the chemical potentials of quarks are assumed to be equal to fi . When g — > , 
only the term A remains. The pressure p ~ A T 4 , with the relations in (45), 
coincides with expression (43). 
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In the particular case of N c — 3 , and taking into account that /i — > as 
T — > oo , we have 

8tt 2 / 21 Ar 
45 I 32 T 

M = -i (l + . 

2 / 1 \ 3 / 2 

Comparing the pressure p = AT A with the Stphan-Boltzmann limit from the 
end of section 5, we get 

- ■ 1 + a 2 g 2 + a 3 g 3 



PSB 

where 



^2 ^3 

B2 = v a » = :v 

Thus, as # — > , the pressure tends to the Stephan-Boltzmann form. 
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Figure Captions 



Fig.l 

The probability of the quark-gluon plasma channel as a function of tempera- 
ture 6 = k B T in MeV and of the relative baryon density n,B/n 0B . 

Fig.2 

The n -meson channel probability on the temperature-baryon density plane. 
Fig.3 

The summary probability of the r] -, p -, and u -meson channels. 
Fig.4 

The nucleon channel probability. 
Fig.5 

The six-quark channel probability. 
Fig.6 

The channel probability of six-quarks in the Bose-condensed state. 
Fig.7 

The quark-gluon plasma channel probability at zero temperature as a function 
of the relative baryon density. 

Fig.8 

The nucleon channel probability at zero temperature vs. the relative baryon 
density. 

Fig.9 

The six-quark channel probability at zero temperature vs. the relative baryon 
density. 

Fig.10 

The pressure (in units of J 4 ) of the multichannel model on the temperature- 
baryon density plane. 

Fig. 11 

The energy density (in units of J 4 ) of the multichannel model. 
Fig.12 

The pressure-to-energy density ratio of the multichannel model. 
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Fig.13 

The specific heat (in units of J 3 ) for the multichannel model. 
Fig.14 

The reduced specific heat of the multichannel model. 
Fig.15 

The compression modulus (in units of J 4 ) of the multichannel model. 
Fig.16 

The glueball channel probability as a function of temperature. 
Fig.17 

Comparison of the glueball and gluon channel probabilities. 
Fig.18 

Reduced specific heat for the gluon-glueball mixture. 
Fig.19 

Relative energy density for the gluon-glueball model without the correcting 
term (solid line), as compared with the lattice Monte Carlo data [42]. 

Fig.20 

Relative energy density for the gluon-glueball model without the correcting 
term (solid line) compared with the lattice numerical simulations [43]. 

Fig.21 

Relative energy density and pressure for the gluon-glueball model satisfying 
the principle of statistical correctness (solid line) as compared with the lattice 
numerical data [43]. 
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Table Caption 



The parameters of the particles included into the numerical investigation of 
thermodynamics of the multichannel model. 
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